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We reconsider the strength of the electroweak phase transition (EWPT) in the inert doublet dark 
matter model, using a quantitatively accurate form for the one-loop finite temperature effective 
potential, taking into account relevant particle physics and dark matter constraints, focusing on a 
standard model Higgs mass near 126 GeV, and doing a full scan of the space of otherwise uncon- 
strained couplings. We find that there is a significant (although fine-tuned) space of parameters for 
achieving an EWPT sufficiently strong for baryogenesis while satisfying the XENONIOO constraints 
from direct detection and not exceeding the correct thermal relic density. We predict that the dark 
matter mass should be in the range 60 — 67 GeV, and we discuss possible LHC signatures of the 
charged and CP-odd Higgs bosons, including a ~ 10% enhancement of the ft — > 77 branching ratio. 



1. INTRODUCTION 

Models of scalar dark matter (DM) H can have inter- 
esting connections to Higgs boson {h) physics because of 
the dimension 4 operator One obvious conse- 

quence is the possibility of the invisible decay channel 
h — >■ HH if the dark matter is sufficiently light. An- 
other is that such a coupling can allow the electroweak 
phase transition (EWPT) to become first order, and po- 
tentially strong enough to be interesting for baryogenesis 
[1] . There has been considerable interest in the interplay 
between dark matter and the electroweak phase transi- 
tion in recent years [2]- [10]. 

The Inert Doublet Model (IDM) is a widely studied 
setting for scalar dark matter [11-14] that can have rich 
phenomenological consequences [15]- [26]. Recently its ca- 
pacity for giving a strong EWPT was considered by rcf. 
[27]. That work found a rather large allowed region of 
parameter space where the EWPT could be strong and 
other constraints satisfied, including the correct thermal 
relic DM density. However, it employed a simplified ver- 
sion of the finite-temperature effective potential, keeping 
only terms up to 0{m/T)^ in the high-temperature ex- 
pansion. On the other hand, a very quantitative treat- 
ment of the effective potential for two-Higgs doublet 
models was recently undertaken in ref . [28] . Our purpose 
in this paper is to reexamine the strength of the EWPT 
in the IDM using this more accurate potential. Moreover 
we search the full parameter space of the model using 
Monte Carlo methods, rather than a restricted subspace 
using a grid search as was done in [27]. We also focus on 
values of the standard model-like Higgs boson mass near 
126 GeV, the value favored by recent LHC data [29]- [32]. 

In this way we are able to extend the results of [27], 
confirming that there exists a significant (though finely 



tuned) region of parameter space in the IDM where 
the strength of the EWPT is sufficiently enhanced for 
electroweak baryogenesis while satisfying other necessary 
constraints. The paper is structured as follows. We re- 
view the definition of the model and collider mass con- 
straints in section 2, our methodology for defining the 
effective potential and scanning the parameter space in 
section 3, and we present the results of the Monte Carlo 
search in section 4. Prospects for testing the model at 
colliders are discussed in section 5, and we give conclu- 
sions in section 6. 



2. THE INERT DOUBLET MODEL 

The Inert Doublet Model is the extension of the stan- 
dard model (SM) by an additional Higgs doublet S with 
the discrete Z2 symmetry S — ^ —S, which naturally leads 
to a stable dark matter candidate in in one of the compo- 
nents of S [12, 14, 33]. Since an unbroken Z2 symmetry 
forbids Yukawa couplings involving S, the second doublet 
interacts with the SM fields only through its couplings to 
the SM Higgs doublet and the gauge bosons. The scalar 
potential of the IDM is given by 
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4 V 2 
Ai(iJt'i/,)(5tJ5,) + X2iH^'H,)iS^^S,) 
[X^H^'H^^S.Sj + h.c] + Xs{S^'S,f (1) 



We assume that S does not acquire a vacuum expection 
value (VEV), so as to keep the Z2 symmetry unbroken. 
The tree-level scalar mass eigenvalues are 
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= ml + ^{Xi + X2 + 2\3)v'^ 
= ml + ^{Xi + X2 - 2X3^ 
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where m^, is the SM-hke Higgs mass, niH is the 

mass of CP-even (odd) component of the inert doublet, 
and m± is the mass of the charged Higgs. Without loss 
of generality, we can take A3 < so that mn < mA and 
therefore H is the dark matter particle. The case A3 > 
just corresponds to renaming H ^ A. We further restrict 
A2 + 2A3 < so that niH < m± to avoid the charged state 
being dark matter. 



2.1. Collider Mass Bounds 



We implemented the relevant phenomenological and 
consistency constraints as in [28], namely precision elcc- 
troweak observables (EWPO), collider mass bounds as 
described in section 2.1, vacuum stability, and the ab- 
sence of Landau poles below 2 TeV (an arbitrary cutoff, 
but sufhcient for considering the model to be a valid ef- 
fective theory up to reasonably high energies). 

For the present study, we add to the above criteria the 
requirement of the correct thermal relic density of dark 
matter. The relic abundance of a dark matter particle H 
is given by [37, 38] 



Precision measurement of the Z boson decay width 
at LEP I forbids the Z boson decay channel Z — > HA, 
which requires that mn + ruA > rnz- In addition, LEP 
II constraints roughly rule out the triangular region [34] 

iTiH < 80 GeV, niA < 100 GeV, - > 8 GeV 

We take the lower bound on the charged scalar mass 
m± > 90 GeV [35]. Following the recent LHC exclusion 
of SM-like Higgs masses in the region 127-600 GeV [29]- 
[32], we restrict m/i to the window 115 — 130 GeV, with 
special attention to the currently favored value m/j = 126 
GeV. 



METHODOLOGY 



3 X 10~27 cm^s-i 



0.1123 ±0.0035 (4) 



Depending on the DM mass itidm, different annihilation 
channels contribute to the thermally averaged annihila- 
tion cross section. Here we consider itidm < rnw so that 
the only relevant annihilation channels are those which 
are mediated by the SM-like Higgs boson into final state 
// pairs (excluding the top quark). Specifically, we con- 
sider the dark matter mass window 45 — 80 GeV (to be 
justified by the results below), and the annihilation cross 
section 
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In this work wc employ the Landau-gauge one-loop 
finite-temperature effective potential similar to that de- 
scribed in ref. [28], which considered the most general 
two-Higgs doublet potential. It has the zero-temperature 
one-loop corrections and counterterms to insure that 
tree-level mass and VEV relations are preserved, and in- 
cludes contributions from the scalars, vectors, Goldstone 
bosons, and the top quark. It further implements resum- 
mation of thermal masses. 

As in [28], we search the full parameter space of the 
model using a Markov chain Monte Carlo (MCMC). 
Models are chosen in such a way as to favor those with 
large values of {vc/Tc), the ratio of the Higgs VEV to 
the critical temperature, which is the figure of merit for 
a strong electroweak phase transition, for the purposes of 
electroweak baryogenesis (see ref. [36] for a review). In 
addition, we favor models with small values of Xdm, the 
effective coupling of the DM to the Higgs boson. 



Ada/ — (Ai + A2 + 2A3) 



(3) 



since this is required by the XENONIOO direct detection 
constraint (see below). Therefore we bias the MCMC us- 
ing the combination a = (vc/Tc)/ Xdm , which is designed 
to produce chains of models such that the probability 
distribution dP/da, treating the chain as a statistical en- 
semble, goes like a. We took the random step size for 
each of the free parameters of the potential to be ^ 10% 
of their starting values, determined by a seed model that 
satisfied the constraints enumerated next. 



where toda/ = = 0.003 GeV is the decay width 

of the Higgs boson (at nih = 126 GeV), Xdm is given by 
(3), and the sum is over all kinematically accessible SM 
fermions, thus dominated by bb pairs in the final state. 

In addition to the WMAP constraints on relic density, 
there is also a strict limit on the spin-independent dark 
matter-nucleon cross section coming from direct detec- 
tion experiments, notably XENONIOO [39]. The relevant 
cross section in the present model is given by [12] 
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where /i = mnruDM / (jnn + 'mDM) is the DM-nuclcon re- 
duced mass. The Higgs-nucleon coupling / is subject 
to hadronic uncertainties that have been discussed in 
refs. [40, 41], and more recently in [42]; in particular 
the quark matrix element ct^at upon which / depends 
is poorly determined. Many authors take / ^ 0.35; for 
example DarkSUSY [43] uses / = 0.38, while ref. [44] 
finds / = 0.35 and a recent estimate based on lattice 
gauge theory [45] obtains / = 0.32 [45]. These do not 
reflect the full range of possible values , which ref. [42] 
puts at / = 0.26 — 0.63, corresponding to a factor of 6 un- 
certainty in the direct detection cross section. We adopt 
the median value / ~ 0.35 for definiteness, but one might 
reasonably invoke / = 0.26 to weaken the effect of direct 
detection limits (by a factor of 1.8 in the cross section) 
and thus the degree of fine-tuning of model parameters 
that we will find below. 
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Because the same Feynman diagram is responsible for 
both processes (5) and (6), the XENONIOO constraint 
restricts \dm to be smalL Thus to get a large enough 
annihilation cross section (5), we must be somewhat close 
to the resonance condition moM — 'm-h/'^, and the model 
thus requires some moderate tuning. We will quantify 
this below. The XENONIOO 90% c.l. limit is cts/ < 8 x 
10~^^ cm^ in the DM mass range of interest, assuming 
that mh = 126 GeV. 

One further condition we impose is that invisible de- 
cays of the SM Higgs boson h — >■ HH do not dominate 
its width, where the invisible contribution is given by 



r. — 

^ inv 



647rm/j 



1 



(7) 



with SM Higgs VEV V = 246 GeV. Ref. [46] indicates 
that a constraint at the level of 40% on the branching 
ratio for such decays should be attainable in the rela- 
tively near future from LHC data. Although not yet 
established, it seems unlikely that the invisible decays 
dominate the width of the Higgs if the excess events seen 
at LHC are really due to the Higgs, so we provisionally 
impose the 40% constraint. Below we will show that this 
requirement restricts the range of allowed toda/, but only 
slightly more than the combination of relic density and 
direct detection constraints. 

We also explore a generalization of the above scenario, 
in which the H boson could make a subdominant con- 
tribution rtn to the total dark matter density. This will 
occur if Xdm is larger, for a given value of uidm, than 
what is required to satisfy (4). The constraint on asi 
from direct detection is correspondingly weakened since 
the rate goes like finer si ^ asi/{<Jv). The factors of 
Xdm cancel out and thus the combined constraints be- 
come independent of Xdai so long as it is large enough 
to sufficiently suppress the relic density. 



4. MONTE CARLO RESULTS 

We initially performed a MCMC scan of the full model 
parameter space (varying Ai, A2, A3, As, ml) as described 
above, without imposing any constraints upon the relic 
density or direct detection cross section. We find in this 
way many models that satisfy the sphaleron constraint 
Vc/Tc > 1 on the ratio of the Higgs VEV to the critical 
temperature. These results are illustrated in fig. 1 which 
plots the DM relic density versus SM-like Higgs mass 
rrih for the models with a strong phase transition (black 
points). It can be seen that there are many such examples 
within the mass window = 115 — 130 GeV, as well as 
within the 3cr allowed range flli^ G [0.085,0.139] for the 
relic density as determined by Seven- Year Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) observations [38], 
and indicated by the shaded horizontal band. A frac- 
tion of the points also satisfy the XENONIOO constraint; 
these are denoted by blue crosses. A small population 
can be found that simultaneously satisfy both constraints 
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FIG. 1; Scatter plot of dark matter relic density flh^ versus 
Higgs mass from Monte Carlo, for models with a strong 
first order electroweak phase transition {vc/Tc > 1). Dense 
points (black) correspond to models that do not necessarily 
satisfy the XENONIOO constraint, while crosses (blue) indi- 
cate models that do. The shaded band corresponds to relic 
density Qh^ £ [0.085, 0.139] observed by WMAP at 3cr. 
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FIG. 2: Scatter plot of \dm versus moM for models with 
strong EWPT, correct relic density (dark points), and ruh ~ 
126 GeV. The 90% c.l. upper bound on Ada/ from XENONIOO 
is shown by the slanted line. The light shaded points de- 
note models whose relic density is subdominant, Q,Hh^ < 
0.085, but which still satisfy the correspondingly relaxed 
XENONIOO limit. The other curves indicate the upper limit 
on Xdm from requiring that the branching ratio for the invis- 
ible decay h — >■ HH not exceed 10% or 40%, respectively. 



and which have rrih near 126 GeV. We use these as seeds 
for a more focused MCMC search in which only mod- 
els with the correct relic density and small enough direct 
detection cross section are admitted into the chains. 

Highlighting the effect of the direct detection con- 
straint, figure 2 shows the scatter plot of MCMC models 
having Vc/Tc > 1 and correct relic density, in the plane of 
Xdm versus thdm for the case of = 126 GeV which 
is suggested by the recent results from the ATLAS and 
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FIG. 3: Monte Carlo distributions of model parameters and 
derived quantities satisfying all constraints (including Q.Hh? £ 
[0.085,0.139]), with mh fixed at 126 GeV. Light-shaded re- 
gions indicate the proportion of frequency contributed by 
models with Vc/Tc < 1, while dark corresponds to Vc/Tc > 1. 
A is the energy scale of the Landau pole for each model. 
Masses, A and the critical temperature Tc are in GeV units. 
The spin-independent DM-nucleon scattering cross section 
asi is in units of cm^ . 



220 



240 



260 



280 



(GeV) 



300 



320 



FIG. 4: Scatter plot of charged Higss mass m± versus GP-odd 
Higgs mass ttia, for models satisfying all constraints. Dark- 
shaded upper region (crosses): models with with Vc/Tc > 1; 
lighter-shaded middle region (squares): models with subdom- 
inant DM component ilnh^ < 0.085 but Vc/Tc > 1; lightest- 
shaded lower region (dots): models with Vc/Tc < 1. The 
diagonal line shows where m± = tua- 



CMS experiments. The XENONIOO upper limit on Ada/ 
is plotted as the slanting line. The peculiar V-shape of 
the allowed region is due to the need for being close to 
resonance of the virtual Higgs boson in the s-channcl 
to get sufficiently strong annihilation for the relic abun- 
dance. Larger values of Xdm do not require to be 
as close to m/j/2. The different density of points above 
and below the Xenon constraint is due to using a rough, 
preliminary MCMC to find the former points, and the 
more focused search to find the latter. 

We also show in fig. 2 the upper limits on Xdm from 
requiring that the branching ratio for the invisible de- 
cays BR{h HH) not exceed 40% or 10% respectively. 
These are futuristic requirements, but ref. [46] estimates 
that the 40% limit will be attainable with 20 fb~^ of in- 
tegrated luminosity at the LHC. When combined with 
the allowed region in our plot this constraint leads to the 
lower bound mjjM ^ 60 GeV. This is only slightly more 
stringent than the bound due to direct detection on the 
left arm of the "V". The interior region of the "V" is 
populated by models for which H makes a subdoniinant 
contribution to the relic density. In this sample, the con- 
straint on the invisible width of the Higgs can be more 
important than the XENONIOO constraint, over a wider 
range of H masses. 

A summary of the favored values of the Lagrangian pa- 
rameters, mass spectra, and EWPT and DM attributes, 
from the refined MCMC search implementing all con- 
straints (assuming H accounts for all of the DM and 
not just a subdominant component), is given by the his- 
togram plots of fig. 3, where the light-shaded portions 
of the bars indicate the fraction of models in the chain 
having Vc/T^ < 1 for the given parameter value, while 
the dark shaded part shows the proportion contributed 



by models with I'c/Tc > 1. The chain contains 14400 
models, of which 2400 have Vc/Tc > 1, despite the bi- 
asing of the MCMC toward large values of Vc/Tc- This 
shows that obtaining a strong first order EWPT is not 
altogether easy. 

Another striking feature is that Ai is always large for 
the cases with Vc/Tc > 1, in the range 2.6 — 3. One might 
expect such large couplings to cause a breakdown of per- 
turbation theory. Nevertheless we have run the renor- 
malization group equations up to find the Landau pole 
energy scale A in each case, keeping only those models 
with A > 2 TeV. The values of A for our accepted mod- 
els arc in the range 2 — 5 TcV; thus some additional new 
physics must come into play at these higher energies. 

A consequence of the large values of Ai and IA2I is that 
there must be fine-tuning between Ai and the combina- 
tion A2 + 2A3 at the level of (1.4 ± 0.5)% (the minimum 
level of tuning we find is 2.5% among the models in our 
chain with Vc/Tc > 1) in order to make Xdm sufficiently 
small. As we noted above, there is an additional tuning 
between 2m dm and mh, which we find to be at the level 
of (3.7 ± 1.4)% in the Vc/Tc > 1 subsamplc. The fine- 
tuning problem for Xdm is slightly ameliorated in the 
scan of models having Vlnh? < 0.085, for which we find 
that A2 + 2A3 must cancel Ai at the level of 3.8 ± 2.7%, 
a factor of 3 improvement. 

5. IMPLICATIONS FOR LHC 

The IDM presents a challenge for discovery at collider 
experiments. From fig. 3, we see that the masses of the 
CP-odd and charged Higgs bosons A and are pre- 
dicted to be in a relatively narrow window 260 — 320 
GeV. We find that there is no strong correlation between 
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niA or m± and m]j]\i, but there is a noticeable corre- 
lation (due to the EWPO constraint) between m± and 
niA, shown in fig. 4. The allowed region consistent with 
Vc/Tc > 1 and H providing the dominant DM component 
is shown by the dark crosses in the upper part of the re- 
gion of interest. This gets extended to somewhat smaller 
masses for the case of models where H is a subdominant 
DM component. 

Because the inert doublet does not couple to quarks, 
its production cross section is suppressed relative to that 
of the SM Higgs, proceeding mainly by qq — > {A,H), 
{H^,H^), {H^,A) or {H^,H) through a virtual Z or 
W in the s channel. Ref. [19] argues that a promising 
discovery signal is through the subsequent decay of the 
heavy Higgs boson(s) to produce lepton pairs and missing 
energy. However the models shown to be most favorable 
for discovery with 100 fb~^ of integrated luminosity arc 
those with |to^ — mojvfl ~ 40 — 80 GeV or muM ^ 
40 GeV, which do not include the ones we predict in 
the present analysis. A similar conclusion holds for the 
trilepton discovery channel [23]. 

Ref. [47] more generally considers monojet events in- 
duced by dimension 6 operators of the form qqxxl ^ for 
Dirac DM x coupling to quarks, producing a jet from 
initial state radiation of a gluon from a quark, as well as 
missing energy from the DM. The IDM generates such 
operators with a virtual Z boson at one loop, through 
the ZZHH coupling, leading to a dimension 6 operator 
of order [G p /^Q''^'^)q$q\H\^ . Even though the effective 
description is not valid at LHC energies, it should give 
an upper bound on the strength of this virtual process. 
Due to the initial state radiation, one of the quarks is 
off-shell, so $ will be of order the momentum of the jet 
rather than the quark mass. Even so, we find that the 
q$q\H\^ operator produces an amplitude similar to that 
of <i<ixxl ^ with A ~ 7 TcV, which is far beyond the 
current sensitivity of LHC experiments, A > 700 GeV, 
determined by ref. [47]. The mismatch arises because of 
the loop suppression factor in our model. The monojet 
constraint is more stringent for operators generated by 
tree-level exchange. For the IDM, the most important 
such process is qq — > HA via s-channel Z exchange, fol- 
lowed by A — >■ HZ. Hadronic decays Z — > qqH could 
then appear as monojets plus missing energy, if the Z is 
sufficiently boosted so that only one jet is resolved. It 
could be interesting to undertake a study of such pro- 
cesses tailored to the IDM. 

Another possibility for testing the IDM is by the effect 
of the charged Higgses on the /i — > 77 branching ratio 
[21, 26]. The rate for /i — > 77 is given by 

where Asm = —6.52 for = 126 GeV and 



with T = {mh/2m±)'^ . The approximation Ah+ = — 1/3 
holds in the limit m± ^ rrih and Aiu^ = 2m j_ which 
are satisfied in the models favored by our analysis. 
Therefore the contribution interferes constructively 
with that of the SM, and results in a increase close to 
10% in the /i — > 77 partial width for all the models 
that we consider. However there is a larger effect on the 
branching ratio for /i — >■ 77 due to the invisible decay 
channel h — >■ HH, which increases the total decay width 
relative to that in the standard model. Thus the change 
in the branching ratio BR(/i — > 77) may be dominated 
by the dilution due to h — > HH in the region where 
niH < nih/'i. If on the other hand mn > mt/'i-, we 
have the definite prediction that BR(/i — > 77) is close to 
110% of its standard model value. Such an enhancement 
is in the direction of the upward fluctuation in the value 
currently observed by CMS [30]. 



6. CONCLUSIONS 

We have quantitatively reconsidered the impact of an 
inert Higgs doublet on the strength of the clcctroweak 
phase transition, generally confirming the result of ref. 
[27] that it is relatively easy to find models satisfying 
all constraints and giving Vc/Tc > 1, but we differ in the 
details. Whereas ref. [27] finds allowed masses m± = itia 
in the range 280 — 370 GeV, we find a more restricted 
range 260 — 320 GeV, correlated with our limit Ai < 3 in 
contrast to theirs, Ai < 4. Considering the approximate 
nature of the finite-temperature effective potential used 
in [27], the results seem to be in reasonable agreement. 

However, we have pointed out some fine-tunings 
needed to make the scenario work: moM must be within 
^ 4 GeV of m/i/2 to get a strong enough HH — > bb anni- 
hilation cross section for the observed relic density, while 
the large value of Ai needed to get Vc/Tc > 1 is canceled 
typically at the (1-4)% level (depending upon whether 
H is the dominant component of the total DM density) 
by A2 4- 2A3 so that the DM-nucleon coupling is small 
enough to satisfy the XENONIOO constraint. If the IDM 
dark matter plus strong EWPT scenario should be borne 
out by experiments, it will be mysterious why these two 
seemingly unlikely coincidences should exist. 

We have considered the prospects for testing the sce- 
nario at the LHC. There are several possible signatures: 
invisible decays of the SM-like Higgs into dark matter 
could be inferred if rriDM ^ 61 GeV. A 10% increase in 
the partial width for h — >■ 77 is a firm prediction. Current 
analyses of missing energy plus monojets or dileptons do 
not seem able to rule out the model in the near future, but 
we suggest for further study that the process qq — > HA 
followed by A — >■ ZH and hadronic decays of the Z (if 
sufficiently boosted) could give monojet-like events that 
might constitute a more promising signal. 
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